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Abstract 

There  is  no  simple  way  to  calculate  the  perimeter  of  an  ellipse.  In  this  expository  paper 
we  review  and  compare  four  methods  of  evaluating  this  perimeter.  These  methods  are 
known  by  the  names  Maclaurin,  Gauss-Kummer,  Cayley,  and  Euler.  Maclaurin’s  method 
is  derived  in  detail,  while  the  other  three  are  simply  described.  A short  introduction  to 
hypergeometric  functions  is  included,  since  these  functions  are  useful  for  finding  the 
elliptic  perimeter.  We  compare  the  usefulness  of  these  four  calculations,  An  appendix 
gives  MICROSOFT  EXCEL®  programs  for  computing  the  perimeter. 
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1.  Introduction 

The  area  of  the  ellipse 


is  given  by  the  formula  A = nab.  We  recognize  this  as  a simple  generalization  of  the  formula 
for  the  area  of  a circle  of  radius  a given  by  A = na2.  What  can  we  say  about  the  perimeter 
of  this  ellipse?  The  circle  of  radius  a has  circumference  C = In  a,  so  is  the  perimeter  of  the 
ellipse  something  simple  like  P = 2n(a  + b)/ 2,  where  we  have  used  the  average  of  a and  b in 
place  of  the  radius  of  the  circle?  The  answer  is  no. 

Unfortunately,  there  is  no  simple  way  to  express  the  perimeter  of  the  ellipse  in  terms  of 
elementary  functions  of  a and  b.  To  express  this  perimeter  we  need  to  expand  our  tool  kit 
of  functions  beyond  the  trigonometric,  exponential,  and  logarithmic  functions  studied  in  the 
calculus.  We  shall  require  ‘special  functions’  such  as  the  gamma  function,  the  elliptic  functions, 
and  the  hypergeometric  functions,  just  to  name  a few.  In  this  paper  we  will  take  a quick  look 
at  ‘complete  elliptic  integrals’  and  ‘hypergeometric  functions’  in  our  quest  to  calculate  the 
perimeter  of  the  ellipse. 

In  Sections  2 and  3 we  will  discuss  in  detail  one  method  of  calculating  the  perimeter  known 
as  the  Maclaurin  expansion  (see  [5]) 


P = 2 na 


f /l\2fc2  / 1 • 3\2&4  / 1 • 3 • 5\2/:6 
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where  k2  = 1 — b2/a2.  This  result  will  also  be  expressed  in  terms  of  elliptic  integrals 
and  hypergeometric  functions  in  Section  4.  In  Sections  5 and  6 we  take  a quick  look  at 
some  properties  of  hypergeometric  functions,  and  in  Section  7 we  introduce  three  additional 
formulas  for  finding  the  perimeter  of  an  ellipse  without  giving  their  derivation.  In  Section  8 we 
use  the  results  of  computer  calculations  to  compare  the  four  formulas  discussed  in  the  paper. 
Appendix  A gives  MICROSOFT  EXCEL  programs  for  computing  the  perimeter. 


2.  The  perimeter  and  elliptic  integrals 

The  differential  arc  length  for  a curve  given  by  parametric  equations  x = x (6)  and  y = y (#) 


is 


ds 


^y+f-Y 

doj  \d0  J 


dO. 


The  ellipse  given  by  the  parametric  equations  x — a cos  0 and  y = b sin  0 has  differential  arc 

length  

d s = yj (— a sin  0)2  + ( b cos  0)2  d 0. 

Therefore,  the  perimeter  of  the  ellipse  is  given  by  the  integral 

ott/2 


= 4/  4a 

JO 


2 sin2  0 + b2  cos2  0 d#, 


in  which  we  have  quadrupled  the  arc  length  found  in  the  first  quadrant.  Replacing  sin2  0 by 
1 — cos2  0 we  get 


rn/z  / 

P = 4 1 y/a2(  1 — cos2  0)  + b2  cos2  0 d 0 

Jo 

rir/2  

= 41  yj a2  — ( a 2 — b2)  cos2  0 d 0. 

Jo 


If  we  let 


we  can  write  (1)  as 


2 b2 

k2  = l 

az 


= 4aT/24i 

Jo 


k2  cos2  0 dO. 


(1) 

(2) 

(3) 


Since  cos#  = sin(7r/2  — #),  for  any  function  f(x)  we  can  write 


pti/2  r 

/ / (cos#)  dO  = 

Jo  Jo 


71/2 


7 r 


/(  sin{  - - # ) )d#. 


With  the  substitution  t = 7r/2  — # we  convert  this  to 


rrc/2  pit/2 

/ /(cos#)d#=  / /(sinf)dL 

Jo  Jo 


Therefore,  we  can  replace  the  cosine  in  (3)  by  a sine  to  obtain 

rn/2 

-k2  sin2 Odd. 


p —4a  I"*  4~\ 

Jo 
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The  reason  for  this  conversion  is  to  express  the  perimeter  in  terms  of  the  named  function 

/^/2  / — 

E(k)=  \j\-k2  sin2  0 dd,  (4) 

Jo 

which  is  called  a complete  elliptic  integral  of  the  second  kind.  Thus  we  can  now  write 

P = 4 aE(k). 

The  complete  elliptic  integral  of  the  first  kind  defined  by  the  integral 

r*/2  dd 


F{k) 


=/ 


y/\  — k2  sin2  0 

occurs  in  the  analysis  of  the  motion  of  the  pendulum  with  large  amplitude  [7,  pp.  331-335]. 

3.  Maclaurin  series  expansion  for  the  perimeter 

Unfortunately  we  cannot  express  the  elliptic  integral  (4)  in  terms  of  elementary  functions. 
To  assist  in  the  integration  of  (4),  we  begin  by  expanding  the  square  root  in  the  integrand  using 
the  binomial  theorem.  We  have 


yfT^X  = (1  ~X)1/2 


(5) 


o°  /lv 

=S0)(-1,V- 


7=0 

in  which  the  binomial  coefficient 

l'p\  P P - 1 P - 2 p-j  + l 


,7/12  3 j 

is  valid  for  j = 1, 2,  3, ... , while  (q)  = 1 • Thus,  for  y = 1, 2, 3, ....  we  have 


i I _ 1 I _ 9 I — 3 

o o 1 2 ^ 2 J 


1 


■7  + 1 


_i  i_l  2-1  3-1 
2 / 2 f 2 2 

1 2 


7-1-5 


We  can  rewrite  this  as 

(— iy 

Returning  to  (5)  we  can  now  write 
\/\—x 


C)-H 


3 4 j 

3 5 2y  — 3 


4 6 8 


2j 


= 1 1 3 5 22-3,, 

2^  2 4 6 8 


7=1 


V 


The  complete  elliptic  integral  (4)  can  be  written  as 


E{k)  = 


13  5 2j  - 3 


2 4 6 8 


2j 


k2j  sin2/ 0 | d6>. 
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Integrating  term  by  term  we  have 


E(k)  = -~Y] 

2 “ 2 4 6 8 

j= 1 


1 1 3 5 2j  -3 

V 


& 


r*/2 

i / sin 

Jo 


in2'  0 d(9. 


(6) 


We  will  make  use  of  the  known  definite  integral  (available  from  software  such  as  MATH- 
EMATICA®  or  standard  mathematical  tables) 

r/2  • 2ioAa  1 • 3 • 5 • • • (2j  — 1)  n 
J0  2 • 4 • 6 • • • 2 7 2 

to  rewrite  (6)  as 


1 4 6 8 

7=1 

2./  - 3 
2 j 

\ /I  3 5 7 
J \2  4 6 8 

27  - 1 
27 

V2-' 

7T 

[i-fY!---.. 

1 A\2  4 6 8 

7=1  x 

27  - 1 

1 \{ 

13  5 7 

27  - 1 

2 

2; 

2j  -l)\2  4 6 8 

2y 

(7) 


This  last  relation  is  the  Maclaurin  series  expansion  for  the  complete  elliptic  integral 


1 

1 

_ f 

1 • 3 \ 2 /c4 

(l  -3-5 \2*6 

“ 

2/  1 V 

2^4/  T" 

\2 • 4 • 6/  5 

4.  Hypergeometric  functions  and  the  perimeter 

We  can  also  express  the  perimeter  of  the  ellipse  in  terms  of  hypergeometric  functions  which 
are  defined  by 

F(a,p;y-,x)  = f^(fJf)jxj,  \x\  < l, 

U {Y)jJ ! 

where  we  have  used  the  Pochhammer  symbol 

(a)  j = a (a  + l)(a  + 2)  • • • (o'  + j - 1). 

Rearranging  (7) 

3 5 2 7 — 3 \ / 1 3 5 7 


Eik)  = | 


^/_1  1 3 5 2y-3W 

p0\  2 4 6 8 2 j A 


2 4 6 8 


27  - 1 

2j 


k2J 


we  rewrite  the  first  factor  in  the  sum  in  (8)  as 


1135  2 j - 3 —2  1 — 2 2 — 2 3 — 2 j-l-l 


2 4 6 8 


2j 


(-A 


(8) 


(9) 
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The  second  factor  in  the  sum  in  (8)  can  also  be  expressed  in  terms  of  Pochhammer  symbols  as 


5.  Hypergeometric  function  identities 

In  this  section  we  summarize  some  facts  concerning  the  important  hypergeometric  functions 
without  giving  their  derivations. 

The  hypergeometric  function  y = F (a,  y \ z)  is  the  solution  of  the  differential  equation 

z(  1 - z)y"  + [y  - (a  + fi  - 1 )z]y'  - afiy  = 0. 

The  following  elementary  functions  can  be  expressed  in  terms  of  hypergeometric  functions: 


1357  2j-\_\\  + \\+2  2+7  — I 


(10) 


Replacing  the  factors  in  (8)  by  (9)  and  (10)  we  can  now  write 


which  can  be  expressed  in  terms  of  the  hypergeometric  function  as 


(11) 


so  that  the  perimeter  of  the  ellipse  is 


(1  +x)p  = F(—p,  1;  1;  -x). 


(1  -x2)1/2  = + 


ln(l  +x)  = xF(  1,  1;  2;  — x), 


We  have  shown  that  the  complete  elliptic  integral  of  the  second  kind  can  be  expressed  in  terms 
of  the  hypergeometric  function  (11),  and  we  can  also  express  the  complete  elliptic  integral  of 
the  first  kind  as 
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The  following  are  known  transformations  for  the  hypergeometric  function: 

f(u,  ($■  2fi;  ^^2)  = (1  + D^F^a,  or  - 0 + 1;  0 + 1;  z2),  (12) 

F(«,  2f}-,  2z)  = (1  — z)-f(|.  a-^,„  + I;  (,3) 

A list  of  notable  mathematicians  who  have  contributed  to  the  theory  of  the  hypergeometric 
function  include  Colin  Maclaurin  (1698-1746),  Leonhard  Euler  (1707-1783),  John  Landen 
(1719-1790),  James  Ivory  (1765-1842),  Carl  Friedrich  Gauss  (1777-1855),  Ernst  Eduard 
Kummer  (1810-1893),  and  Arthur  Cayley  (1821-1895). 


6.  Generalized  hypergeometric  functions 

In  this  section  we  briefly  introduce  the  generalized  hypergeometric  function  defined  by  the 
series 

F ( a\,a2,  ...,am  \ _ {a\)k(ai)k  • • • ( flmh  xk 

mtn\  bx,b2,...,bn  X ) ~ L (bi)k(jb2)k  ■ ■ ■ (bn)k  k\' 

If  m < n the  series  converges  for  all  x,  if  m = n + 1 the  series  converges  for  \x\  < 1,  and  if 
m > n + 1 then  the  series  diverges.  Another  notation  in  common  use  is 


mFn 


( a\,a2,  ...,0m 

V b\>b2,-.-,bn 


xj  = mFn(a\,a2,  b\,b2, bn\x). 

The  ordinary  hypergeometric  function  considered  in  previous  sections  can  thus  be  written  as 
F(a , b\  c\  x ) = 2F\  ^ = 2F\  ( a , b\  c\  x). 

All  of  the  formulas  listed  in  Section  5 without  derivation  can  be  found  in  [1]. 


7.  Three  additional  formulas  for  the  perimeter 


In  this  section  we  list  three  additional  formulas  for  the  perimeter  of  the  ellipse  without  giving 
their  derivation.  Two  of  the  three  can  be  expressed  conveniently  in  terms  of  the  hypergeometric 
function. 

Our  first  expression  is  the  Gauss-Kummer  relation 


P = jt(a  + b) 


1-3-5 
2 • 4 • 6 • 8 


2 

/*8  + -.- 


where 

a — b 

h = -. 

a + b 

We  note  that  the  Gauss-Kummer  relation  is  derived  from  the  hypergeometric  series  [2],  [4] 


P = n{a  + b)F^-1-, -X--\-h2y  (14) 

The  Gauss-Kummer  relation  (14)  is  obtained  by  setting  a = — f and  z = h = 

(a-b)/(a  + b)  in  (12). 
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Our  second  formula  was  derived  by  Cayley  [3,  pp.  46-55]  for  the  case  when  a lb  is  large  by 
substituting  k — b/a  in  the  governing  differential  equation.  Cayley’s  series  is 


Our  third  formula  was  given  by  Euler  and  is  represented  by  the  hypergeometric  relation 

/•=W2(«’  + *)f(I-I;1;  (^)  ).  (16) 

Euler’s  formula  (16)  is  obtained  by  setting  a = — ft  = and  2z  = k2  = 1 — b2  /a2  in  (13). 


8.  Comparing  the  four  formulas  for  the  perimeter 

We  now  use  a computer  to  compare  the  following  four  formulas  for  the  perimeter  of  the 
ellipse: 

/ 1 1 9\  9 b2 

Maclaurin  P = 2anF\  — , 1;  kz  ),  kl  — 1 

V 2 2 / a 2 


Table  1 . 


a 

b 

Perimeter 

Best 

formula 

Gauss- 

Kummer 

Cayley 

Euler 

Maclaurin 

1 

0.9999 

6.282  871  152 

G-K 

2 

76027 

2 

4 

1 

0.99 

6.251808  848 

G-K 

3 

1 192 

4 

9 

1 

0.9 

5.973160433 

G-K 

6 

136 

7 

18 

1 

0.8 

5.672333  578 

G-K 

7 

68 

10 

29 

1 

0.7 

5.382  368  981 

G-K 

9 

44 

14 

42 

1 

0.6 

5.105  399773 

G-K 

10 

31 

20 

62 

1 

0.5 

4.84422411 

G-K 

13 

24 

28 

93 

1 

0.4 

4.602622519 

G-K 

16 

18 

42 

147 

1 

0.3 

4.385  91007 

Cayley 

21 

14 

73 

260 

1 

0.2 

4.202008  908 

Cayley 

30 

11 

153 

562 

1 

0.1 

4.063  97418 

Cayley 

56 

8 

549 

2027 

1 

0.05 

4.019425  619 

Cayley 

102 

6 

1 941 

7131 

1 

0.04 

4.013143  313 

Cayley 

124 

6 

2907 

10646 

1 

0.03 

4.007  90945 

Cayley 

159 

6 

4 880 

17  791 

1 

0.02 

4.003  83916 

Cayley 

225 

5 

10074 

36454 

1 

0.01 

4.001  098  33 

Cayley 

405 

4 

34183 

121721 

1 

0.005 

4.000309233 

Cayley 

723 

4 

112  844 

393  112 

1 

0.004 

4.000205  049 

Cayley 

870 

4 

164534 

568  193 

1 

0.003 

4.000120518 

Cayley 

1 100 

4 

265  848 

906354 

1 

0.002 

4.00056807 

Cayley 

1527 

3 

515  404 

1719  206 

1 

0.001 

4.000155  88 

Cayley 

2 643 

3 

1 520644 

4816476 

1 

0.000  1 

4.00000202 

Cayley 

13  943 

3 

20939  398 

43144616 

1 

0.000001 

4.00000000 

Cayley 

80  386 

2 

- 

- 
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Gauss-Kummer 


Cayley 


a h 

•3/4  2 1 

— In--  — - — 


• 4 V k 1-2  3-4 


l2  • 32  • 5 
+ 22  - 42  -6 

/ 4 2 2 1 

x ( In 


Euler 


P =71 


In  Table  1 we  list  our  results  for  a standard  ellipse  with  semi-major  axis  a = 1 and  semi-minor 
axis  b ranging  from  1 to  0.  The  perimeter  is  listed  in  the  third  column.  In  the  last  four  columns 
we  give  the  number  of  terms  needed  in  the  series  to  compute  the  perimeter  numerically  to  eight 
decimal  places. 


A user-defined  function  in  MICROSOFT  EXCEL  is  easy  to  develop.  To  open  the  spreadsheet, 
press  Alt-Fl  1.  You  will  see  Microsoft  Visual  Basic  Book  1 . Click  ‘Insert’  and  select  ‘Module’. 
You  get  the  Module  1 screen.  Type  in  the  functions  in  visual  basic  as  given  below.  The  functions 
will  be  available  in  the  spreadsheet.  The  three  functions  developed  are  hgs  (a,b,  c,  z) , 
gauskum(a, b) , cayley  (a,  b) , and  gkc  (a,  b) . The  first  one,  hgs  (a,  b) , calculates  the  hy- 
pergeometric series  F(a,  b\  c\  z),  the  other  three  evaluate  the  ellipse  perimeter  (gauskum  ( a , b) 
using  the  Gauss-Kummer  series  (14),  cayley (a, b)  using  the  Cayley  formula  (15),  and 
gkc  (a,  b)  using  the  Gauss-Kummer  or  Cayley  series  so  that  the  perimeter  is  evaluated  using 
the  least  number  of  terms  in  the  series). 

Function  hgs  (a,  b,  c,  z) 
t = a*b/c*z 
s = 1 + t 
n = 1 
Do 

a = a + 1 

b = b + 1 

c = c + 1 

n = n + 1 

t = t*a/c*b/n*z 
si  = s 
s = s + t 

Loop  While  si  <>  s 

hgs  = s 

End  Function 

Function  gkc (a,  b) 

a = Abs ( a ) : b = Abs ( b ) 

If  a < b Then 

c=a:b=a:c=b 


Appendix  A.  MICROSOFT  EXCEL  implementation 
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End  If 

If  b = 0 Then 
gkc  = 4 * a 
Exit  Function 
End  If 

If  b / a > 0.25  Then 
gkc  = gauskum(a,  b) 

Else 

gkc  = cayley (a,  b) 

End  If 

End  Function 

Function  gauskum(a,  b) 
h = (a  - b)  / (a  + b) 
t = 0.5  * h 
p = 1 + t * t 
u = 1 
Do 

t = t * u / (u  + 3)  * h 
pl  - P 

p = p + t * t 
u = u + 2 

Loop  While  pl  <>  p 

gauskum  = Application . Pi ( ) * (a  + b)  * p 
End  Function 

Function  cayley (a,  b) 
x = b / a 

t2  = Log (4  / x)  - 0.5 
x = x * x 
tl  = 0 . 5 * x 
p = 1 + tl  * t2 
al  = 1:  a2  = 2 
si  = al  / a2 
s2  = 1 / al  / a2 
Do 

tl  = tl  * si 
t2  = t2  - s2 
al  — al  + 2 
a2  = a2  + 2 
si  = al  / a2 
s2  = 1 / al  / a2 
tl  = tl  * x * si 
t2  = t2  - s2 
pl  = p 

p = p + tl  * t2 
Loop  While  pl  <>  p 
cayley  = 4 * a * p 
End  Function 


The  calculation  terminates  when  the  addition  of  a term  does  not  change  the  value  for  given 
machine  precision.  The  precision  is  about  le-12. 
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Michon  [6]  discussed  the  Gauss-Kummer  and  Cayley  series;  the  article  provides  the 
QBASIC  function  gk  (h)  using  a combination  of  Gauss-Kummer  and  Cayley  formulas. 
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